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Photons have proven to be excellent carriers of quantum information, and play essential roles in 
numerous quantum information processing (QIP) applications jJ-Qj]. In particular, heralded single¬ 
photon sources [5!] via spontaneous parametric-down conversion (SPDC) have been a key technology 
for demonstrating small-scale QIP, yet their low generation efficiency is a critical limitation for further 
scaling up optical QIP technology. In order to efficiently overcome the probabilistic nature of SPDC, 
here we demonstrate time multiplexing flHjjl] for up to 30 time slots of a periodically pumped heralded 
single-photon source, using a switchable low-loss optical storage cavity. We observe a maximum 
single-photon probability of 38.6±0.4% in periodic output time windows, corresponding to ^6 times 
enhancement over a non-multiplexed source, but with no increase in the contribution of unwanted 
multi-photon events. Combi ning this time-multiplexing technique with a heralded source producing 
pure single-photon states | Ill 4l3l| should enable larger scale optical QIP systems than ever realized. 

Realizing and scaling up optical QIP systems requires on-demand preparation of quantum states of light such as 
single-photon and definite multi-photon states. Single atoms, ions, and solid-state single-emitter sources such as color- 
centers in diamond and semiconductor quantum dots can generate true single-photon states, and especially solid-state 
systems have great potential for integration. However, while high indistinguishability [TdJ, [l5j] and somewhat high (50- 
80%) collection efficiency M have been reported in different single-emitter systems, achieving both simultaneously 
remains a challenge. Moreover, most single-emitter sources generate single photons with narrow bandwidths (5-100 
MHz) that may be unsuitable for pursuing high-speed applications (e.g., a 5-MHz bandwidth source can produce 
non-overlapping single-photon gaussian wavepackets—of duration 90 ns—at a maximum rate of ~10 MHz). 

SPDC is another approach that has been conventionally and widely used for generating entangled photon pairs, 
multi-photon entangled states, and small-scale quantum algorithms (5J. A photon pair generated via SPDC can also 
be used for generating a “heralded” single-photon state; detecting one of the photons “heralds” the presence of the 
other. Moreover, current technology has realized photons with very high coupling efficiency into single-mode optical 
fibers EDI and controlled two-photon spectral shape mm- However, the photon-pair generation process via 
SPDC is probabilistic: one cannot obtain a photon pair or a heralded single photon on demand. Unfortunately, 
indefinitely increasing the mean number of photon pairs per pump pulse p also increases the likelihood of unwanted 
fc-photon pairs (~p fc ) by higher-order processes. 

To overcome the probabilistic nature of photon-pair generation, time-multiplexing techniques were first proposed 
and demonstrated by Pittman, Jacobs, and Franson Q in 2002. The method was since extended and theoretically 
analyzed [Ml (a related approach uses spatial multiplexing [l9l - l23| . but this is much more resource intensive— a 
source analogous to ours would need 30 crystals and detectors). The basic idea is shown in Fig. 1. (a). A laser pulse 
train with period r pumps a y*- 2 -* nonlinear crystal, generating photon pairs (i.e., signal and idler photons) in one 
or more time slots. Each signal photon is sent to a single-photon detector (SPD) whose firing heralds in which time 
slot the corresponding idler photon is present. By using a switchable storage cavity with the matched cycle length r, 
any one of the time slots heralded to contain an idler photon can be multiplexed onto a single output time window 
(see Fig. 1 (b)). Thus, the multiplexed single-photon probability Pm{ 1 ) during the output time window is increased 
according to the number of pump pulses (time slots) N used for one cycle of the multiplexing. Moreover, if N is large, 
the probability of generating unwanted multiple pairs in a given time slot can be made arbitrarily small, because the 
total pump energy through the multiplexing cycle is distributed over the N time slots, and the ratio of the single- and 
multi-photon probability is as low as the one for a single (non-multiplexed) heralded single-photon source. Hence, 
assuming losses can be kept low, the multiplexed heralded single-photon source can work as a pseudo-on-demand 
single-photon source (i.e., Pm{ 1) — > 1) by using sufficiently low pump pulse energy (p < 0.1) and large cycle number 
N. A detailed mathematical calculation of Pm( 1), accounting for loss and inefficiencies of SPDs, PC, and the SPDC 
source is discussed in urn and in the Supplementary Information. 

A schematic diagram of our experimental setup is shown in Fig. 2 (a). We implemented the setup with a low-loss 
heralded single-photon source and photon storage cavity (see Methods). However, in practice, many repeated passes 
through the storage cavity attenuate the stored idler photons; for the case that multiple time slots within N are 
heralded, photons stored later—closer to the output time window—experience less loss. In order to select only the 
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FIG. 1: (a) Simplified schematic diagram of a time multiplexed heralded single-photon source, (b) Timing diagram of pump, 
signal, and idler photons in the time multiplexing. In our system, when the idler photon enters the storage cavity via a 
polarizing beam-splitter (PBS), a Pockels cell (PC) in the cavity is fired, rotating the photon polarization by 90°, thereby 
storing the photon in the cavity loop until its polarization is switched back. The single-pass delay of the cavity is set to match 
r so that the time slots where idler photons may be present are temporally and spatially overlapped with each other in the 
cavity. The idler photon is then released from the cavity in the desired output time slot by a second switching of the PC. 
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FIG. 2: Schematic diagram of experimental setup. PBS1: Brewster-angled polarizing beam splitter. PBS2: cube polarizing 
beam splitter. NPBS: non-polarizing beam splitter. SPD: single-photon detector. PC: Pockels cell. IF: interference filter 
(AA = 20 nm). HWP: half-wave plate. QWP: quater-wave plate. FPGA: field-programmable gate array. Dashed red line 
shows the optical path when the optical delay is shortened to AtDL ~200 ns. 


last-produced, and therefore, lowest-loss photon, we introduced a low-loss Herriott-cell-type optical delay line (25| 
(between the SPDC crystal and the switchable storage cavity) whose delay A tuL is long enough to contain idler 
photons from all N of the pump pulses. The long optical delay makes it possible to determine the time slot of the 
last-detected signal photon, and to only switch the corresponding idler photon into the storage cavity. For our largest 
delay, A toL is ~400 ns, sufficiently longer than the ~120-ns electronics latency (from detecting signal photons to 
the PC driver) to hold idler photons from up to N ~30 (30 x r = 250 ns for r = 8.33 ns in our experiment). Using 
ffippable mirrors, we can also reduce A tjjL to ~200 ns, which is still longer than the electronics latency, but not 
enough to store the latest heralded time slot for larger N. This allows us to compare the effect of storing the photons 
in the latest and earliest heralded time slots. 

We studied time multiplexing with two different pump powers: A high pump power (such that p = 0.35) is used so 
that the heralding signal rate H is nearly saturated, i.e., H is close to the repetition rate of the multiplexed source 
R = 50 kHz for N = 30 (see blue circles and red triangles in Fig. 3 (a)). With a low pump power (p = 0.07), 
H increases almost linearly with N up to 30 (see green squares and purple diamonds in Fig. 3 (a)). Note that R 
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FIG. 3: Experimental results, (a) Heralding rate H, (b) single-photon probability Pm( 1), and (c) the second-order correlation 
function </ 2 l(0) versus number of multiplexed time slots TV. Blue circles: p = 0.35, A tDL ~400 ns. Red triangles: p = 0.35, 
A toL ~200 ns. Green squares: p = 0.07, A tDL ~400 ns. Purple diamonds: p = 0.07, A tDL ~200 ns. Error bars are 
estimated by poissonian photon counting statistics. For (a) and (b), solid and dashed curves show theoretical predictions 
for the measurement with p = 0.35 and 0.07, taking into account the measured experimental parameters ( p,T c ,Ti,r/Ds ) and 
using Eq. (3) in Supplementary Information. Solid and dashed lines in (c) indicate the values of g < ' 2 \ 0) measured without 
multiplexing (TV = 1) for p = 0.35 and 0.07. 


for this proof-of-principle multiplexed source is limited by the duty cycle of the PC drivers. Figure 3 (b) shows our 
experimental results for single-photon-in-periodic-output-time-slot Pm(1) versus TV. Pm( 1) is estimated by 


Pm{ 1) 


Si + S2 

RriDi 


( 1 ) 


where S 1 and S 2 are the single-photon count rates measured by SPD1 and SPD2 at the output time windows, 
and rjui is the system detection efficiency for the multiplexed photons (see Methods). For all the measurement 
in our experiment, we clearly observed an enhancement in Pm( 1) above the single-cycle probability: We observed 
Pm( 1) = 38.6 ± 0.4% for p = 0.35 and TV = 30 (see blue dots in Fig. 3 (b)), corresponding ~ 6 times enhancement 
over the non-multiplexed heralded photons (for TV = 1, Pm( 1) = 6.8±0.4%). The enhancement factor was even larger 
for p = 0.07 (~16, see green squares in Fig. 3 (b)), because H increases as TV up to 30. To our knowledge, these 
single-photon probability and enhancement factors are superior to all previous demonstrations mini of multiplexed 
heralded single-photons (that achieved Pm( 1) < 1% with enhancement factor of 1.6-4). For the results with the 200-ns 
delay line (red triangles and purple diamonds in Fig. 3 (b)), as expected Pm( 1) starts to decrease after TV ~15, due 
to increased loss for photons stored in the earlier heralded time slots. Our results are in excellent agreement with our 
theoretical prediction (solid and dashed curves in Fig. 3(b)), except the observed single-photon probabilities with the 
200-ns delay line for TV > 20. The difference may be caused by slight spatial-mode mismatching between the output 
of the delay line and the switchable storage cavity—a 1% mismatch would be sufficient to make slight beam clipping 
at the PC in the storage cavity for TV > 20 and to account for the observed discrepancy. 

We characterized the ratio of single- and multi-photon probabilities via the second-order correlation function at 
zero time delay between the signals from SPD1 and SPD2 |2q : 


<? (2) ( 0 ) 
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where C 12 is the coincidence count rate between SPD1 and SPD2. The observed g^ (0) < 1 (see Fig. 3 (c)) is expected 
for a non-classical light source. Moreover, while the single-photon probability is greatly enhanced by multiplexing, the 
observed g^ (0) is nearly constant versus TV for the same p: the contribution of the two-photon probability relative 
to the single-photon one is approximately independent of the time multiplexing. As expected (solid and dashed lines 
in Fig. 3 (c)), g i ' 2 \ 0) is higher for p = 0.35, because of the higher likelihood of multi-pair emissions compared to the 
case for p = 0.07. Nevertheless, our results are still better than what one could achieve even from an ideal heralded 
SPDC source (e.g., with 100%-efhcient trigger detection and perfect coupling into lossless optics, see Table I), which 
predicts Pm{ 1) = 36.8% (25.0%) for a completely mixed (pure) single-photon state according to poissonian (thermal) 
photon-number statistics [27| . 

There are several straightforward improvements for the time-multiplexed source. First, we observed a low purity 
(0.05) of the heralded idler-photon state (using the method proposed in H3), because of strong spectral entanglement. 
Since many QIP applications utilize multi-photon interference effects, which require pure states |28(, our SPDC source 
must be changed to one that heralds single photons in a pure quantum state, i.e., without spectral entanglement 
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TABLE I: Comparision between non-multiplexed sources and multiplexed sources. The ideal non-multiplexed sources assume 
photon-pair sources with perfect heralding efficiency and 100%-efficient detectors with photon-number-resolving functionality. 
For the time-multiplexed source with efficient components, the values are predicted by feasible parameters (p = 0.1, T c = 99%, 
Ti = 95%, r\os = 90%) using recent state-of-the-art technologies and Eq. (SI) in Supplementary Information. 


Pm( 1) 

fl (2) (0) 


Ideal non-multiplexed source Multiplexed source 

Mixed state Pure state This work With efficient components 


36.8 % 
0 


25.0 % 
0 


38.6 ± 0.4 % 
0.479 ± 0.028 


> 80% 
< 0.05 


between the signal and idler photons. Sources that herald intrinsically pure single photons without narrowband 
filtering, have been reported !ll{l3[; combined with temporal multiplexing, these have the potential for both high 
single-photon probability and purity. (Note that the probability of higher-order emissions can actually be higher for 
pure heralded single-photon sources than for ones generating mixed states, because of their different photon-number 
statistics, i.e., thermal for a single-mode SPDC source versus poissonian for multimode states (27]; multiplexing is 
consequently even more important to suppress unwanted multiple-photon events). 

Beyond preparing pure photons, it is important to maintain this purity throughout the time-multiplexing process. In 
particular, group velocity dispersion (GVD) in the storage cavity needs to be negligibly small, lest the cycle-dependent 
GVD disturb the spectral-temporal mode purity. Thus, more careful dispersion management in the storage cavity is 
required for heralded photons with shorter pulse durations (< 100 fs). Finally, SPDs with higher detection efficiencies 
would better detect unwanted multiple photons in the signal mode, improving the likelihood of heralding true single¬ 
photon states. Recent state-of-the-art detector technologies, such as visible-light photon counters [29| and arrays of 
superconducting nanowire detectors [ffij have achieved > 90% detection efficiencies and timing resolutions below 100 
ps. 

With feasible values for collection (95%) and detector efficiency (90%), and delay-line (1%) and switchable cavity 
losses (1%), we anticipate single-photon probability up to 80%, with a </ 2 )(0) as low as 0.05 (see Table I). Despite that 
our current source rate is only 50 kHz, such performance would already be advantageous for building up larger photon 
numbers; for example, the generation rate of an 8-photon state would be 0.8 s x 50 kHz~8000 per second by using the 
multiplexing technique, while recent 8-photon experiments [fj with non-multiplexed SPDC sources generated only ~1 
per second. 

In conclusion, we experimentally demonstrated a time-multiplexed heralded single-photon source. Our experiment 
is implemented by a highly efficient SPDC source, low-loss fixed optical delay line, and adjustable storage cavity 
with fast polarization switches. We observed large enhancement in single-photon probabilities by multiplexing up to 
30 time slots of the heralded single-photon source. Despite the great enhancement in the single-photon probability, 
the observed second-order correlation function of the output photons showed that the ratio of single-photon and 
unwanted two-photon probabilities for the multiplexed source was as low as for a non-multiplexed one. We anticipate 
that incorporating a bright and pure heralded single-photon source and state-of-the-art photon-number-resolving 
detectors will enable a periodic and near-deterministic single-photon source that will be a critical component for 
larger scale optical QIP applications. 


Methods 

We used a frequency-tripled mode-locked Nd:YAG laser (A = 355 nrn, At ~5 ps, t = 8.33 ns) to pump a bismuth 
barium borate (BiBO) crystal that generates non-collinear photon pairs centered at 710 nm via Type-I phase-matching. 
Each photon produced by the SPDC was passed through an interference filter (IF, A A = 20 nm), and then coupled 
into a single-mode fiber (SMF), with a coupling efficiency above 70% 0- We used a cascade of four SPDs, each a 
silicon avalanche photodiode (Si-APD), with detection efficiency of ~70%. The detector cascade allows us to reduce 
the effect of each SPD’s saturation and to resolve the approximate number of photon pairs produced in a time slot. By 
measuring the ratio of single- and coincidence-count rates from the SPDC source, we estimate that the total system 
detection efficiency of the signal photons rj^s = 41.8 ± 0.4%. A field-programmable gate array (FPGA) module 
processes input signals from the SPDs, and triggers PCs at heralded and output time slots. A Herriott-cell-type 
optical delay line 1251 is made of two concave spherical mirrors separated by ~1 m. The idler photon enters the cell 
through a hole in one mirror, bounces 54 times between the two mirrors, and escapes from the cell through a hole 
in the other mirror. We used one way (round trip) of the Herriott cell to implement ~200 ns (^400 ns), with 

transmission of 89.3% (84.6%). For an adjustable photon storage cavity, we used a custom Brewster-angled polarizing 
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beamsplitter (PBS1) and PC comprising a pair of rubidium titanyl phosphate (RTP) crystals. The single-pass cavity 
transmission (T c = 97.0%) is limited by imperfect reflection of mirrors and PBS1, and the PC’s imperfect transmission 
and polarization switching. A PC and PBS2 after the cavity are used as an optical shutter, which opens in an output 
time window of the multiplexing only if at least one heralding signal is detected within N multiplexed time slots, 
and drops unheralded photons. The total transmission of other optics for the idler mode Ti is 92.4%. In order to 
investigate the multiplexed single-photon probability using Eq. 0> we estimated the system detection efficiency for 
the multiplexed photon r)pn = 56.8%, accounting for the measured transmission of the non-polarizing 50:50 beam 
splitter (NPBS, 89.9%), fiber coupling efficiency and transmission (90.2%), and SPD detection efficiency (70%, as 
listed in manufacturer specification sheet). 


Acknowledgement 

This work was supported by NSF Grant No. PHY 12-12439, US Army ARO DURIP grant No. W911NF-12-1-056, 
and the MURI Center for Photonic Quantum Information Systems (ARO/ARDA Program DAAD19-03-1-0199). 


[1] Gisin, N., Ribordy, G., Tittel, W. & Zbinden, H. Quantum cryptography. Rev. Mod. Phys. 74, 145195 (2002). 

[2] O’Brien, J. L. Optical quantum computing. Science 318, 1567 (2007). 

[3] Childs, A. M., Gosset, D. & Webb, Z. Universal computation by multiparticle quantum walk. Science 339, 791-794 (2013). 

[4] Pan, J.-W. et al. Multiphoton entanglement and interferometry. Rev. Mod. Phys 84, 777-838 (2012). 

[5] Hong, C. & Mandel, L. Experimental realization of a localized one-photon state. Phys. Rev. Lett. 56, 58-60 (1986). 

[61 Pittman, T., Jacobs, B. & Franson, J. Single photons on pseudodemand from stored parametric down-conversion. Phys. 
Rev. A 66, 042303 (2002). 

[7] Jeffrey, E., Peters, N. A. & Kwiat, P. G. Towards a periodic deterministic source of arbitrary single-photon states. New J. 
Phys. 6, 100 (2004). 

[8] McCusker, K. & Kwiat, P. G. Efficient Optical Quantum State Engineering. Phys. Rev. Lett. 103, 163602 (2009). 

[9] Glebov, B. L., Fan, J. & Migdall, A. L. Deterministic generation of single photons via multiplexing repetitive parametric 
downconversions. Appl. Phys. Lett. 103, 031115 (2013). 

[10] Mower, J. & Englund, D. Efficient generation of single and entangled photons on a silicon photonic integrated chip. Phys. 
Rev. A 84, 052326 (2011). 

[11] Mosley, P. J., Lundeen, J. S., Smith, B. J. & Walmsley, I. A. Conditional preparation of single photons using parametric 
downconversion: a recipe for purity. New J. Phys. 10, 093011 (2008). 

[12] Evans, P. G., Bennink, R. S., Grice, W. P., Humble, T. S. & Schaake, J. Bright source of spectrally uncorrelated 
polarization-entangled photons with nearly single-mode emission. Phys. Rev. Lett. 105, 253601 (2010). 

[13] Spring, J. B. et al. On-chip low loss heralded source of pure single photons. Opt. Exp. 21, 13522-13532 (2013). 

[14] Sipahigil, A. et al. Indistinguishable Photons from Separated Silicon-Vacancy Centers in Diamond. Phys. Rev. Lett. 113, 
113602 (2014). 

[15] Wei, Y.-J. et al. Deterministic and robust generation of single photons from a single quantum dot with 99.5% indistin- 
guishability using adiabatic rapid passage. Nano Lett. 14, 6515-6519 (2014). 

[16] Gazzano, O. et al. Bright solid-state sources of indistinguishable single photons. Nature Commun. 4 (2013). 

[17] Christensen, B. G. et al. Detection-loophole-free test of quantum nonlocality, and applications. Phys. Rev. Lett. Ill, 
130406 (2013). 

[18] Pereira, M. C. et al. Demonstrating highly symmetric single-mode, single-photon heralding efficiency in spontaneous 
parametric downconversion. Opt. Lett. 38, 1609 (2013). 

[19] Migdall, A. L., Branning, D. & Castelletto, S. Tailoring single-photon and multiphoton probabilities of a single-photon 
on-demand source. Phys. Rev. A 66, 053805 (2002). 

[20] Ma, X., Zotter, S., Kofler, J., Jennewein, T. & Zeilinger, A. Experimental generation of single photons via active multi¬ 
plexing. Phys. Rev. A 83, 043814 (2011). 

[21] Collins, M. J. et al. Integrated spatial multiplexing of heralded single-photon sources. Nature Commun. 4, 2582 (2013). 

[22] Francis-Jones, R. J. A. & Mosley, P. J. Exploring the limits of multiplexed photon-pair sources for the preparation of pure 
single-photon states. arXiv quant-ph 1409.1394v2 (2014). 

[23] Bonneau, D., Mendoza, G. J., O’Brien, J. L. & Thompson, M. G. Effect of loss on multiplexed single-photon sources. New 
J. Phys. 17, 043057 (2015). 

[24] McCusker, K. T. Efficient quantum optical state engineering and applications. Ph.D. Thesis, Appendix C (2012). 

[25] Robert, C. Simple, stable, and compact multiple-reflection optical cell for very long optical paths. Appl. Opt. 46, 5408-5418 
(2007). 

[26] Loudon, R. The Quantum Theory of Light (Oxford Univ. Press, Oxford, 2000). 

[27] Christ, A., Laiho, K., Eckstein, A., Cassemiro, K. N. & Silberhorn, C. Probing multimode squeezing with correlation 
functions. New J. Phys. 13, 033027 (2011). 



6 


[28] Hong, C., Ou, Z. & Mandel, L. Measurement of subpicosecond time intervals between two photons by interference. Phys. 
Rev. Lett. 59, 2044-2046 (1987). 

[29] Kim, J., McKay, K. S., Kwiat, P. G., Zielnicki, K. & Gansen, E. J. Novel Semiconductor Single-Photon Detectors. 
Experimental Methods in the Physical Sciences 45, 147-183 (2013). 

[30] Marsili, F. et al. Detecting single infrared photons with 93% system efficiency. Nature Photon. 7, 210-214 (2013). 


Supplementary Information 

Here we show theoretical details of a time-multiplexed single-photon source with inefficient optical components. As 
discussed in the main text, for the case that multiple time slots within N are heralded, it is the more efficient to use 
the photon that is created last. The probability of emitting an m-photon state by multiplexing N time slots and 
storing the last-born photons is given by 

N oo oo 

P M {m) =5^(1 - ]T J2 Pc(k 2 )P d (l\k 2 )P e (m\j,k 2 ), (3) 

3=1 ki=l fc 2 =l 

where P c (k), P d (l\k), and P e {m\j,k) are, respectively, the probability of generating k photon pairs in a time slot, 
detecting a signal photon conditioned by k photon pairs created, and emitting m idler photons conditioned by k 
photon pair generations in the j -th time slot. The first part in Eq. © is the probability that a single signal photon 
detection has not occurred in the “last” j — 1 time slots. The second part is the probability that a single signal photon 
is detected at the j-th time slot, and to photons are released at the output time window. For pure heralded photons 
via SPDC generating separable photon pairs, P c (k) follows a thermal distribution [27]: P c (k) = p k /(1 + p) k+1 , where 
p is the mean number of created photon pairs per time slot. However, the observed purity of our heralded idler 
photons is Trp 2 ~0.05, indicating that our SPDC source generates photon pairs into s ~ 1/Trp 2 ~20 Schmidt modes. 
Assuming that each of s Schmidt modes has the identical mean photon numbers per time slot p/s, P c (k ) follows the 
distribution: 


P c (k) = 


(P/s) 


k + s — 1 
k 


(1 +p/s) k+s 

For s>l, and thus for our SPDC source, Eq. © closely follows a Poissonian distribution: 


P c (k) ~ 


e -Ppk 

k\ ' 


The detection probability of a signal photon P d (l\k) is given by 


P d {l\k) = J2 r lDs( l -VDs) k l (^j (Jj 


l-l 


( 4 ) 


( 5 ) 


( 6 ) 


where r/Ds is the total transmission of the signal photons from the SPDC crystal to SPDs, and D is the number 
of SPDs used for a trigger-detector cascade to herald idler photons. Here, we assume that the SPDs are “bucket” 
detectors that discriminate between zero and one-or-more photons, and the detector cascade distributes the signal 
photons to D detectors with an equal probability (1/D). By substituting D —> oo, Eq. © can also describe a 
photon-number-resolving detector: P d (l\k) = r]Ds(l — r]D S ) fc_1 . The idler photon’s transmission to an output time 
slot P e (m\j, k) is given by 


P e (m\j, k) = (Hlf-*) m (l - Ti T?-i) k - m . (7) 

The T c and T) denote the idler photon’s storage efficiency in the switchable storage cavity and other optics transmission 
(including an initial delay line, fiber coupling efficiency, etc.), respectively. For the case without an extra optical delay 
line to determine the last heralded idler photon, Pm (to) for storing the earliest-born photons is given by (7[: 

N oo oo 

P M (m ) =51(1- 55 PcfaWIW' -1 51 Pc(k 2 )Pd(l\k 2 )P e (m\j,k 2 ). 
j =i fci=i fc 2 =i 


( 8 ) 
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The above equation is similar to Eq. ©) except that the first part in Eq. (J5J describes the probability that the single 
signal photon is not detected in the “first” j — 1 time slots. As shown in Fig. 3 (c) in the main text, Eq. © has 
an optimal number of time slots -/V max to maximize Pm^), because increasing N effectively results in larger photon 
storage times, and associated larger loss in the storage cavity. -/V max is obtained by solving <9 Pm (w)/<9iV = 0 in N: 

M _ log(log(l -a)) -log(logr c ) 
max " log T c — log(l — a) ’ 

OO 

a = ^ PMPaillh). (10) 

ki = l 

For our experiment, N max ~14 and ~29 for p = 0.35 and 0.07, respectively. 



